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SUMMARY 

The  laminar  boundary  layer  on  a  circular  cone  at  angle  of  attack 
to  a  supersonic  stream  is  analyzed.  A  solution  for  the  outer  nonviscous 
flow,  as  obtained  by  a  method  of  perturbation  in  angle  of  attack,  is 
available  and  is  presumed  to  govern  the  development  of  -the  boundary 
layer,  subject  to  certain  restrictions  and  corrections; 

In  the  limit  of  vanishing  cone  angle,  the  boundary- layer  equations 
reduce  to  those  for  the  flow  about  a  yawed  infinite  cylinder  in  com¬ 
pressible  flow,  but  may  not  be  solved  in  the  simple  manner  appropriate 
to  incompressible  flow. 

For  small  angle  of  attack,  the  boundary-layer  quantities  are  con¬ 
sidered  to  differ  only  slightly  from  the  corresponding  quantities  in 
the  known  flow  at  zero  angle  of  attack.  The  system  of  equations  result¬ 
ing  from  this  assumption  are  solved  to  yield  the  effect  of  small  angle 
of  attack  on  the  velocity  profiles,  skin  friction,  boundary-layer  thick¬ 
ness,  and  the  forces  and  moment  on  the  cone,  or,  more  precisely,  the 
rates  of  change  of  these  quantities  with  angle  of  attack,  evaluated  at 
zero  angle  of  attack. 

All  boundary-layer  quantities  show  parabolic  similarity  in  merid¬ 
ional  planes.  Circumferentially,  the  boundary- layer  thickness  associated 
with  the  profile  of  the  circumferential  velocity  component  is  constant, 
and  the  profile  itself  exhibits  no  tendency  to  separate;  whereas  the 
profile  of  the  meridional  component  undergoes  changes  caused  by  angle 
of  attack  which  show  the  expected  effects  of  secondary  flow  in  the 
boundary  layer. 


INTRODUCTION 

The  flow  around  bodies  of  revolution  at  positive  angles  of  attack 
is  characterized  by  a  draining  of  the  boundary  layer  from  the  underside 
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to  the  top,  causing  the  formation  of  "lobes"  of  low-energy  air  near  the 
top  of  the  "body.  For  sufficiently  large  angles  of  attack,  separation  of 
the  cross  flow  may  result  in  a  transversely  shed  Karman  vortex  street. 
Under  these  circumstances,  large  viscous  contributions  to  lift,  drag, 
and  pitching  moment  on  the  body  may  be  expected.  Of  course,  viscous 
lift,  drag,  and  moment  are  expected  at  any  angle  of  attack,  whether  or 
not  cross-flow  separation  occurs. 

The  present  investigation,  conducted  at  the  MCA  Lewis  laboratory, 
is  concerned  with  a  special  case,  namely,  the  laminar  boundary  layer 
which  develops  on  a  right  circular  cone  at  a  small  angle  of  attack  to 
the  incident  supersonic  flow.  Solutions  of  the  laminar  boundary- layer 
equations  are  sought  in  order  to  describe  the  viscous  effects  due  to 
angle  of  attack,  and  particularly  those  contributing  to  lift,  drag,  and 
pitching  moment. 

As  is  customary,  it  will  be  assumed  that  the  nonviscous  solution 
for  the  supersonic  flow  about  a  cone  at  angle  of  attack  is  valid,  except 
in  a  thin  viscous  layer  enveloping  the  cone  surface.  The  nonviscous.  flow 
is  conical  in  the  Busemann  sense;  that  is,  physical  quantities  are  con¬ 
stant  along  rays  from  the  apex.  According  to  the  analysis  of  Stone 
(reference  l) ,  the  flow  quantities  may  be  considered  to  be  the  sums  of 
the  corresponding  quantities  of  zero  angle  of  attack  and  small  increments 
linear  in  angle  of  attack.  This  theory  forms  the  basis  for  the  extensive 
tables  prepared  at  Massachusetts  Institute  of  Technology  under  the  super¬ 
vision  of  Kopal  (references  2  and  3).  The  tables  of  references  2  and  3 
will  be  used  in  the  present  report.  In  reference  4,  Ferri  criticizes 
the  Stone-Kopal  treatment  of  the  entropy  gradients  near  the  body  and 
proposes  a  correction  to  the  tables  of  reference  3.  The  applicability 
of  this  correction  will  be  discussed  in  a  subsequent  section. 

The  boundary  layer  itself  is  to  be  treated  in  the  manner  of  refer¬ 
ence  5,  wherein  it  is  shown  that  "conical"  bodies  in  supersonic  flow 
have  boundary  layers  which  develop  parabolically  along  generators  of  the 
surface .  This  boundary- layer  development  is  a  consequence  of  the  fact 
that  the  outer  nonviscous  flows  over  such  bodies  have  vanishing  velocity 
and  pressure  gradients  along  rays  from  the  apex. 

The  special  case  of  a  right  circular  cone  at  zero  angle  of  attack 
has  been  solved  by  Hantsche  and  Wendt  in  reference  6  and  by  Mangier  in 
an  unavailable  report.  These  solutions  give  a  boundary  layer  which 
thickens  parabolically  along  cone  generators,  and  provide  a  formula  for 
viscous  drag  which  is  valid  also  for  small  angle  of  attack,  as  will  be 
shown  subsequently.  '  •  • 

A  comparison  is  made  herein  between  the  present  analysis  and  the 
hypothesis,  which  has  been  discussed  by  Allen  (reference  7)  and  Van  Dyke 
(reference  8),  that  the  viscous  cross  flow  associated  with  the  flight 
of  bodies  of  revolution  at  angle  of  attack  be  regarded  as  essentially 
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that  which  would  occur  in  the  case  of  an  infinite  circular  cylinder  at 
yaw  with  respect  to  the  incident  stream,  a  flow  which  has  been  analyzed 
by  Prandtl,  by  R.  T.  Jones,  and  by  Sears  (reference  9).  The  attractive¬ 
ness  of  this  hypothesis  lies  in  the  fact  that,  for  the  yawed  infinite 
cylinder  in  incompressible  flow,  the  cross  flow  may  be  obtained  from 
two-dimensional  theory  without  regard  to  the  spanwise  component.  In 
reference  5,  however,  it  is  pointed  out  that  the  chordwise  velocity  com¬ 
ponent  does  depend  to  some  degree  on  the  spanwise  component  in  the  case 
of  compressible  flow  about  a  yawed  infinite  cylinder.  Whether  or  not 
this  consideration  is  of  practical  importance  has  not  yet  been  ascertained. 


NONVISCOUS  FLOW  ABOUT  CONE  AT  ANGLE  OF  ATTACK 

In  order  to  treat  the  boundary  layer  on  a  cone  at  angle  of  attack, 
it  is  necessary  first  to  describe  the  nonviscous  flow  about  the  body, 
particularly  the  velocity  components  and  properties  of  state  which  would 
be  predicted  at  the  surface.  Figure  1  shows  the  coordinate  system  and 
other  notation  to  be  used.  The  distance  x  is  measured  along  generators 
of  the  cone;  y,  normal  to  the  cone  surface;  and  <p,  circumferentially 
around  the  cone.  The  velocity  components  u,  v,  and  w  are  in  the 
x-,  y-,  and  <P -direction,  respectively.  The  cone  has  a  semivertex 
angle  of  0  and  is  at  an  angle  of  attack  a.  The  subscript  0  denotes 
conditions  in  the  undisturbed  stream,  and  subscript  1  denotes  eval¬ 
uation  of  the  nonviscous  flow  at  the  cone  surface,  or,  alternatively, 
at  the  outer  edge  of  the  boundary  layer.  The  symbol  notation  used  herein 
appears  in  appendix  A. 

Since  the  present  report" is  concerned  with  small  angles  of  attack, 
and  correspondingly  small  departures  of  the  entire  flow  from  that  occur¬ 
ring  at  zero  angle  of  attack,  it  will  be  convenient  to  refer  all  physical 
quantities  to  this  basic  nonviscous  flow.  Thus,  denoting  by  bars  the 
properties  at  zero  angle  of  attack,  evaluated  at  the  cone  surface,  the 
following  symbols  on  the  left  will  hereinafter  be  identified  with 
the  dimensionless  groups  on  the  right: 

u,  w  *  u/u,  w/u 

p  ~  p/p 

T  -  2CpT/u2 

P  ~  p/pu2 
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Thus,  for  example,  T 
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On  the  surface  of  the  cone,  the  theory  for  small  angle  of  attack 
used  in  reference  3  yields  results  which  may  he  written  in  the  following 
form: 

U]_  =  1  -  a  Aj  cos  <p 


w^  =  a  Ar,  sin  <P 
p/p  =  1  +  a  A3  cos  (p 


> 


pi 


=  1  +  a  Aa  cos  <p 


(2) 


J 


where  the  A!s  correspond  as  follows  to  the  quantities  (evaluated  at 
the  cone  surface)  tabulated  in  reference  3  and  do  not  themselves  depend 
on  a: 


Aq 

-x/u 

a2 

-z/u  -  2x/ (u  sin  <3>) 

a3  ~ 

Tl/P 

a4  - 

?/p 

J 


(3) 


The  barred  (a  =  0)  quantities  are  tabulated  in  reference  2 .  The 
quantities  on  the  right  side  of  equations  (3)  are  in  the  symbol  notation 
of  reference  3. 


The  relations  (2)  provide  the  boundary  conditions  to  be  imposed  on 
the  viscous  flow  at  the  outer  edge  of  the  boundary  layer. 

If  the  definitions  used  in  reference  3  were  to  be  followed  strictly, 
then  Ag  ~  z/u.  The  corresponding  relation  in  set  (3)  is  written  dif¬ 
ferently  in  order  to  compensate  for  a  systematic  error  in  the  tables. 

This  error,  confirmed  by  Professor  Kopal  in  a  letter  to  the  Lewis 
laboratory,  consists  in  the  use  of  the  wrong  sign  for  the  quantity  C 
appearing  in  equation  (34)  of  reference  3.  Only  the  values  of  z  (and 
hence  A2)  are  affected.  Since  C  refers  to  the  deviation  of  entropy 
in  the  flow  due  to  the  cone  thickness,  the  correction  described  is  of 
little  consequence  for  small  cone  vertex  angles,  but  is  substantial  for 
the  larger  vertex  angles. 

In  reference  4,  it  is  pointed  out  that  the  relations  (2)  imply  that 
the  change  in  entropy  at  the  cone  surface  due  to  angle  of  attack  varies 
circumferentially  around  the  cone  as  cos  whereas  physical  reasoning 


* 


w 


m 
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indicates  that  at  the  cone  surface  the  entropy  should  he  constant  except 
for  a  singularity  at  =  jt .  This  difficulty  arises  because  the  method 
of  linearization  in  angle  of  attack  used  in  reference  3  is  improper  near 
the  cone  surface.  In  reference  4  it  is  concluded  that  the  tables  of 
reference  3  are  correct,  except  in  a  "vortical  layer"  of  an  angular 
thickness  of  order  a2  (infinitesimally  thin  in  the  linear  approximation) 
across  "which  there  are  a  sudden  change  in  entropy  of  order  a  and 
corresponding  changes  in  other  flow  quantities,  such  as  velocity  and 
temperature,  but  not  pressure.  The  "vortical  layer"  is  a  transition 
region  in  which  the  form  of  the  entropy  variation  must  change  from  the 
type  cos  <p  away  from  the  surface  to  constant  at  the  cone  surface. 

The  corrections  to  the  tables  of  reference  3  required  to  take 
account  of  this  effect  are  not  used  in  the  present  report;  that  is, 
the  boundary  layer  will  be  considered  as  governed  by  the  flow  external 
to  the  layer  alluded  to  in  reference  4.  The  following  argument  is  pre¬ 
sented  to  justify  this  point  of  view:  The  reason  for  considering  the 
entropy  constant  at  the  surface  is,  briefly,  that  the  fluid  particles 
nearest  the  surface  must  all  have  passed  through  the  conical  bow  shock  wave 
in  the  vicinity  of  <p  =  0  and,  once  past  the  shock,  will  conserve  their 
entropy  at  the  value  appropriate  to  the  strength  of  the  shock  at  <p  =  0. 
Thus,  the  sharp  normal  gradients  near  the  surface  (discussed  in  refer¬ 
ence  4)  may  be  considered  to  exist  by  reason  of  the  conservation  of 
entropy  along  the  streamlines  nearest  the  wall.  If  the  boundary  layer  is 
much  thicker  than  the  "vortical  layer"  (a2),  the  boundary  layer  will 
immediately  entrain  most  of  the  fluid  in  the  vicinity  of  <p  =  0,  which 
would  otherwise  proceed  along  the  surface  to  form  the  base  of  the  vor¬ 
tical  layer.  The  properties  of  the  entrained  particles  will  then  dif¬ 
fuse  through  the  boundary  layer  by  the  action  of  viscosity,  entropy  will 
no  longer  be  conserved  next  to  the  wall,  and  the  constraint  giving  rise 
to  the  vortical  layer  will  be  removed. 

Thus,  it  is  reasonable  to  suppose  that  viscous  diffusion  in  the 
boundary  layer  will  prevent  the  formation  of  the  vortical  layer  described 
in  reference  4,  and  that  the  boundary  layer  will  be  governed  at  its 
outer  edge  by  the  flow  described  in  reference  3  (equations  (2)),  provided 
that  the  boundary  layer  is  much  thicker  than  would  be  the  vortical  layer. 

If  the  reverse  were  true,  the  boundary  layer  would  certainly  be  governed 
by  conditions  more  nearly  appropriate  to  the  base  of  the  vortical  layer. 

The  criterion  for  entrainment  and  diffusion  of  the  vortical  layer  would 
be 

6  >>xot,2 

where  5  is  the  boundary- layer  thickness.  As  a  rough  approximation 
(reference  6), 
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■where  R  Is  Reynolds  number  per  unit  length.  Thus,  the  criterion  is 

Rx«1-4 

Since  the  present  analysis  will  concern  only  small  angles  of  attack,  and, 
strictly  speaking,  will  he  precise  only  in  the  limit  a-^0,  consistency 

requires  that  the  inequality  Rx<<ar4  he  presumed  applicable .  ^  Of 
course,  since  the  vortical  layer  is  associated  with  entropy  variation, 
this  question  is  of  less  importance  for  smaller  cone  vertex  angles . 

LAMINAR  BOUNDARY  LAYER  ON  CONE  AT  ANGLE  OF  ATTACK 
Differential  Boundary-Layer  Equations 
and  Boundary  Conditions 

In  reference  5  it  is  shown  that  the  differential  equations  for  the 
compressible  laminar  boundary  layer  on  a  conical  body  in  supersonic 
flow  may  be  written: 

jf  +  |  Ell£l  g  +  |  g gj  f\\  -  |  gXfX£  +  |  (gx>2  +  2fXXX  =  0  (4a) 

0 

(4b) 

T  +  (fx)2  +  (g\)2  =  T1  +  ui2  +  wi2  (4c) 

p  =  riipT  (4d) 

Equations  (4a)  and  (4b)  are  momentum  equations,  equation  (4c)  is  an 
energy  balance,  and  equation  (4d)  is  the  equation  of  state.  The  func¬ 
tions  f(X,£)  and  g(\,£)  are  related  to  the  tvo -component  vector 
potential  discussed  in  reference  5  and  are  defined  according  to  the 
relations 


w  =  gxJ 


and  in  a  manner  such  as  to  satisfy  the  continuity  equation  identically. 


[f  - 1  ^ e  + !  gJ 


g 


XX 


3  sxsx£ 


2  „  -  2  p 

3  SX  X  “  3 


in 


+  2SXXX  = 
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These  equations  embody  the  following  physical  assumptions,  in 
addition  to  those  arising  from  the  concept  of  a  thin  boundary  layer 
(for  example,  that  pressure  does  not  vary  across  the  boundary  layer)  : 

(^a)  Prandtl  number  of  1  and  constant  specific  heats 

(b)  No  heat  transfer  through  the  body  surface 

(c)  Validity  of  the  temperature -viscosity  relation 

t.C?  (6; 


where  C  is  defined  as  follows,  in  order  to  match  relation  (6a)  to  the 
Sutherland  formula  at  the  cone  surface  (denoted  by  subscript  w) : 


■(f) 


1  +  S/T 
Tw/T  +  S/T 


the  quantity  S  being  taken  equal  to  (216°  R)  2Cp/u2.  This  temperature- 

viscosity  relation  as  applied  to  flow  over  a  flat  plate  is  discussed  by 
Chapman  and  Rubesin  in  reference  10.  Since  the  case  of  Prandtl  number 
of  1  and  no  heat  transfer  is  under  consideration,  Tw  may  be  taken  equal 

to  the  stream  stagnation  temperature  (equation  (4c)). 

The  coordinates  x  and  y  (fig.  l)  have  been  made  dimensionless 

by  referring  them  to  the  length  C  ==.  Thus,  the  normalizing  relation 

pu 


x,y  -  e®,  £S 

cjr  Cii 

should  be  added  to  the  conventions  (l)  previously  established 
The  coordinate  X  has  been  formed  as  follows: 


■^§) 


t  rv 


p  dy  x  2 


Equations  (4)  and  definition  (8)  show  that  similarity  of  the  Blasius  type 
exists  in  meridional  planes.  Thus,  for  example,  boundary- layer  thickness 
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is  proportional  to  xl/2  and  may  be  said  to  increase  parabolically 
with  x.  This  fact  is  a  consequence  of  the  vanishing  pressure  gradient 
along  generators  of  a  conical  body  in  supersonic  flow. 


The  variable  £  is  a  dimensionless  coordinate  which,  in  the  case 
of  a  circular  cone,  may  be  replaced  by  <f>  sin  0,  or  by  <p0  where 
0  =  sin  0.  Equations  (4)  may  thus  be  rewritten: 


g  + 


-sl 

30  S<Pj 


XX 


2 


30 


V2  +  2fxxx  “  0 


(9a) 


r  1  p'(cp)  2  "1  2  2  _  2  p'(<p)  Q  n 

[5  +  30  p  g  +  30  SXX  “  30  SXgX<p  "  3  sXfX  '3?  p  +  2gXXX  “  0 

(9b) 

T  +  (f\)2  +  (g\)2  =  +  u-l2  +  Wj_2  (9c) 

p  =  Ili  pT  (9d) 

The  boundary  conditions  on  the  functions  f(\,<p)  and  g(\,<p) 
are:  At  the  outer  edge  of  the  “boundary  layer ,  the  u  and  w  velocity 
components  should  take  on  the  corresponding  nonviscous  values 

^(•j^)  =  Uqtv)  (10a) 

gx(-><j>)  =  wi(<p)  (1Gb) 

and  at  the  cone  surface,  the  u,  v,  and  w  velocities  should  vanish, 

f\(0,<p)  =  gx(0,<p)  =  o  (10c) 

f(0,<P)  =  g(0,<p)  =  0  (lOd) 

Equation  (lOd)  is  equivalent  to  the  requirement  that  the  normal  velocity 
component  v  vanish  at  the  cone  surface.  This  boundary  condition  is 
discussed  at  some  length  in  reference  5. 


Simplification  of  Equations  in  Special  Limiting  Cases 

The  problem  of  finding  solutions  of  equations  (9)  for  f  and  g 
subject  to  boundary  conditions  (10)  is  rather  difficult.  Accordingly, 
it  is  desirable  to  investigate  the  manner  in  which  the  equations 
simplify  in  certain  special  circumstances . 
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Small  cone  angle ,  large  angle  of  attack.  -  Symbolically,  this  case 
may  be'  described  as  follows: 


®«1 

a 


a  of  order  1 


These  -assumptions  would  he  expected  to  lead  to  the  equations  for  the 
boundary  layer  on  a  yawed  infinite  cylinder  (the  compressible  analog 
of  the  flow  treated  in  reference  9) .  For  the  purpose  of  this  dis¬ 
cussion,  it  is  convenient  to  define  new  variables: 


A  s  yb  *  (I) 


P  <ly 


(6x) 


(U-) 


f  h  t/^fze 

.G  =  g/^30 


Equations  (9)  and  (10)  become 


30F  +  2^*1.  G  +  2G<p]  FM  -  2GAFA(p  +  20  (Ga)2  +  2F ^  =  0  (12a) 


30F  +  G  +  2Gj  GM  -  2GAGA(p  -  20GAFA  -  2  +  20^  =  0 


(12b) 


Fa(-,<P)  =  u-l(<p) 

ga(”^)  =  wi(<p) 

FA(0,q>)  =  GA(0,<p)  =  F(0,<p)  =  G(0,q>)  =  0 


(13) 


For  flow  about  a  cone  at  angle  of  attack,  F  and  its  derivatives  will  be 
presumed  to  be  of  unit  order;  G  and  its  derivatives,  and  the  derivatives 
°f  Pi*  of  order  a  (order  1,  in  this  case) .  When  terms  of  order  0 
are  neglected,  equations  (12)  may  be  written: 


GAFAcp  -  GJpFAA  -  P2p<P'>  gfAA  +  fAAA  (14a) 

^AGA<P  "  G®GAA  =  "  P  (-1-  ”  GGAA)  +  gAAA  (14b) 
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Equations  (14)  are  also  the  equations  for  the  compressible  viscous 
flow  about  a  yawed  infinite  cylinder  and,  for  the  incompressible  case, 
are  precisely  equivalent  to  those  treated  in  reference  9.  Inspection 
of  definition  (ll)  shows  that  similarity  exists  in  terms  of  the  reciprocal 
square  root  of  the  local  radius  (0x)  of  the  cone  cross  section.  For 
an  infinite  cylinder  at  yaw,  this  radius  would  be  a  constant.  For  a 
sufficiently  slender  cone  at  a  sufficiently  large  angle  of  attack 
(0/a«l),  the  hypothesis  discussed  in  reference  7  is  therefore  essen¬ 
tially  correct;  namely,  that  at  any  point  on  the  cone,  the  boundary 
layer  is  that  which  would  appear  for  the  same  angle  of  attack  on  an 
infinite  circular  cylinder  having  a  diameter  equal  to  the  cross- 
sectional  diameter  of  the  cone  at  the  point  under  consideration. 

Equations  (14a)  and  (14b)  are  coupled  through  the  density  p,  which 
appears  in  the  pressure  gradient  term,  and  therefore  G  cannot  be 
obtained  independently  of  F  as  in  the  incompressible  case.  Compress¬ 
ibility  enters  in  this  way  according  to  the  stream  Mach  number.  In 
reference  8,  it  is  suggested  that  the  important  Mach  number  is  that 
based  on  the  cross  velocity. 

This  limiting  case  retains  the  general  feature  of  parabolic  growth 
of  the  boundary  layer  along  cone  generators  despite  the  fact  that,  in 
the  limit,  the  configuration  is  that  of  a.  cylinder  of  infinite  length. 


Small  angle  of  attack,  large  cone  angle, 
specified  as 

—  <<1 
0 


This  situation,  may  be 


0  of  order  1 


presumably  yielding,  in  the  limit,  the  equations  for  the  boundary  layer 
on  a  cone  at  zero  angle  of  attack  to  a  supersonic  stream.  Equations  (2) 
show  that  for  small  angles  of  attack,  ux(<p)  *  1  and  w-^)  and  p'(<P) 

are  of  order  a.  Therefore,  from  equations  (10),  f(X,<p)  will  he  of 
order  1  and  g(X,<p)  will  be  of  order  a.  With  terms  of  order  a  in 
equations  (9a),  (9b),  and  (10)  neglected,  equations  (9a)  and  (10)  become 


«xx  +  2f\xx  =  0 


(15a) 


fx(->*)  =  1 


(15b) 


fx(0,<p)  =  f(0,<p) 


(15c) 


Equation  (9b)  drops  out  because  all  its  terms  are  of  order  a  or  smaller 
In  view  of  the  form  of  equations  (15),  f  is  a  function  only  of  X,  and 
equation  (15a)  may  be  interpreted  as  an  ordinary  differential  equation. 


2286 


MCA  TN  2521 


11 


Equations  (15)  govern  the  boundary  layer  on  a  cone  at  zero  angle  of 
attack  to  a  supersonic  stream  and  are  formally  identical  with  the  equa¬ 
tions  for  plane  flow  over  a  flat  plate ,  the  difference  between  the  cone 
and  flat  plate  cases  appearing  only  through  a  slightly  different  definition 
of  independent  variable.  ( See,  for  example,  reference  6  or  Mangier' s 
work . ) 


PERTURBATION  ANALYSIS  FOR  SMALL  ANGLE  OF  ATTACK 

In  the  remainder  of  this  report,  the  laminar  boundary  layer  on  a 
cone  at  angle  of  attack  to  a  supersonic  stream  will  be  analyzed  as  a 
perturbation  of  the  flow  at  zero  angle  of  attack  (equation  (15));  that 
is,  under  the  conditions 


9  of  order  1 


Formulation  of  Perturbation  Equations 

With  the  assumption  that  the  entire  flow  may  be  represented  as  the 
sum  of  the  basic  flow  fo(X)  satisfying  equations  (15)  and  correction 
terms  proportional  to  angle  of  attack,  f(X,q>)  and  g(X,<p)  may  be 
written 


f  5  f0M  “  a  Aq  cos  +  •  •  •  (16a) 

g  =  a  A2  sin  <p  g1(X)  +  .  .  .  (16b) 

The  quantities  fg,  f^,  g^,  A^_,  and  Ag  are  presumed  of  order  1.  The 

forms  of  the  correction  terms  are  chosen  to  be  consistent  with  the 
relations  (2)  to  be  imposed  as  boundary  conditions  on  the  boundary  layer 
at  its  outer  edge. 


The  pressure  gradient  terms  appearing  in  equations  (9a)  and  (9b) 
may  be  expressed  in  perturbation  form  with  the  aid  of  equations  (2), 
(9c),  and  (9d)  and  definitions  (16):  From  equations  (2),  since  pressure 
is  presumed  not  to  vary  across  the  boundary  layer. 


p'Gp) 

p 


-a  A3  sin  q> 


(17  a) 


From  equations  (2),  (9c),  and  (9d), 
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i  r-i  t  r-i  T1  +  U12  +  V12  ~  (fx)2  "  (g\)2 

p  2y  p  2y  p  (l  +  a  A3  cos  <p) 


(17b) 


Because  T]_  differs  from  T  only  by  a  quantity  of  order  a,  equa¬ 
tions  (16)  may  be  substituted  into  equation  (17b)  to  yield,  for 
small  a. 


P  _  r-i 


p  2 y 


T  +  1  -  (f0’)  +  .  .  . 


which  may  be  multiplied  into  equation  (17a)  to  give 

- p  — -  =  -a  A3  sin  q>  T  +  1  -  (fQ * )2  +  .  .  •  (17c) 

Substituting  equations  (16)  and  (17c)  into  (9a)  and  (9b)  and 
equating  the  sum  of  terms  of  unit  order  in  equation  (9a)  to  zero  give 


Vo"  +  2fom=  0 


The  sum  of  terms  of  order  a  in  equation  (9a)  is  next  equated  to  zero., 
yielding 

+  3f„f  "  +  3f  "f  +  6f _,n  =  0  (19) 


0  °1  0  1 


0  1  1 


The  largest  terms  in  equation  (9h)  are  of  order  a;  the  sum  of  these 
terms  is  equated  to  zero : 


3W  -  2foV  +  ^i’”  -  -2  e4  ^  T 


H- 


<v> 


Substituting  equations  (16)  and  (2)  into  boundary  conditions  (10) 
and  equating  terms  of  like  order  in  a  provide  the  boundary  conditions 
to  vhich  the  differential  equations  (18),  (19),  and  (20)  are  subject: 


V(-)  =  1 


f  '(0)  =  fo(0)  =  0 


f1'(~)  =  1 


^'(0)  =  fx(o)  =  0 


&!_'(•)  =  1 


(0)  =  gx(0)  =  0 


(23) 
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In  view  of  boundary  conditions  (2l)  and  (23),  the  relation 

^5  r-i  t  =  i 

©Ag  2 


(24a) 


must  hold  in  order  that  equation  (20)  he  satisfied  at  X  =  “>.  Because 
of  relations  (l). 


2  _1_ 
^M2 


(24h) 


where  M  is  the  Mach  number  at  the  surface  of  the  cone  when  at  zero 
angle  of  attack.  In  terms  of  quantities  tabulated  in  reference  3, 


2  u^ 

r- 1  i-u2 


(24c) 


With  the  aid  of  equations  (3),  (24b),  and  (24c),  it  may  be  verified  from 
reference  3  that  equation  (24a)  is  indeed  correct.  Equation  (20)  may 
therefore  be  written  more  simply  as  follows : 


3f0g"  -  2f  -g  -  +  6g  •"  =  -2  <1  +  Izi  m2  1  -  (f  ') 


To  summarize,  the  basic  flow  at  zero  angle  of  attack  over  the  cone 
is  determined  by  the  solution  of  the  boimdary- value  problem 

foV  +  2VM  =  0  (18) 

V(~)  =  1  f0'(°)  =  f0(°)  =  0  (21) 

The  meridional  velocity  component  u  is  given  by  the  function  f0’(X), 
according  to  equation  (5) .  When  fg(X)  is  known  and  M  is  specified, 
the  circumferential  velocity  component  due  to  a  small  angle  of  attack 
may  be  determined  from  equations  (5)  and  (16b)  and  from  the  solution 
of  the  problem 


3fng  "  -  2f, 


6gl’M  =  -2  Si  +  Ili  M2  1  -  (f  ') 


(•)  =  1 


gl’(°)  =  Sx(0)  =  0 


Then,  presuming  Tq(X)  and  g^(X)  to  have  been  determined  and  the 
quantity  Ar,/6A^  to  be  specified,  the  increment  in  meridional  velocity 
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component  due  to  a  small  angle  of  attack  may  be  found  from  equations  (5). 
and  (16a)  and  from  the  solution  of  the  problem 


3f0f1!l  +  3^"^  +  6fi 


(  ^2\ 

[ 11 1 

\@A 1/0  1 


fl'(-)  =  1 


^'(0)  =  f±(0)  =  0 


The  perturbation  procedure  used  to  obtain  these  differential  equa¬ 
tions  is,  in  general,  subject  to  an  important  limitation  which  should 
be  examined.  Since  all  velocity  components  vanish  at  the  cone  surface, 
it  is  perhaps  not  obvious  that  it  is  proper,  near  the  surface,  to  make 
the  neglections  in  equations  (9)  necessary  to  arrive  at  equations  (19) 
and  (25).  In  particular,  consider  the  term  ff^  appearing  in  equa¬ 
tion  (9a) .  Substitution  of  equation  ( 16a)  yields 

ff\X  =  f0f0"  -  a  A-l  cos  (P  (f±f0"  +  f0f]_")  +  a 2  A±2  cos2  <p  f^"  +  .  .  . 

(26a) 


In  deriving  equation  (19),  the  term  involving  or  has  been  neglected 
in  comparison  with  the  term  involving  a,  despite  the  fact  that  both 
terms  vanish  as  \“*0.  The  problem  of  proper  neglection  near  the  sur- 
face  becomes  determinate  if  the  last  two  terms  of  expression  (26a) 
are  written  as  follows : 


w 


-a  A±  cos  <p  (fpfo"  +  fofl")  ll  -  a>  Ax  cos  <p 


f  f  ”  +  f  f  " 
10  0  1  , 


Thus,  the  question  is  whether  or  not,  for  a«  1,  the  following  inequality 
holds  as 


a  A-^  cos  <p 


£lV 

Vo"  +  fofi' 


(26b) 


or,  equivalently,  whether  or  not 


f  f  "  +  f  f  " 
10  0  1 


(26c) 


is  always  of  order  1  or  smaller . 
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Boundary  conditions  (2l)  and  (22)  indicate  that  fQ  and  f-j_  may 
be  represented  for  small  X  by  the  following  power  series: 

fo  =  FT  fo"(o)  x2  +jrfo,,,(°)  x3  +  •  •  • 
fi =  fi"(0)  x2  +  wr  fi"(°)  x5  + .  .  . 

If  f0"(0)  /  0  (as  in  the  present  problem),  then,  in  the  limit  X-*0, 

quantity  (26c)  is  equal  to  i  f i"(0)/fo"(0),  which  is  presumably  no 

2 

larger  than  unit  order.  All  the  neglections  made  in  deriving  equa¬ 
tions  (19)  and  (25)  may  be  justified  by  this  sort  of  argument. 

If,  however,  fQn(o)  were  zero  and  f-^"(0)  were  not  zero,  expres¬ 
sion  (26c)  would  become  infinite  as  X-*0,  and  inequality  (26b)  would  not 
hold.  The  quantities  f0"(0)  and  ^"(o)  are  related  to  the  slopes 

at  the  wall  of  the  basic  and  incremental  velocity  profiles  (equa¬ 
tion  (5)),  and  hence  are  proportional  to  the  corresponding  contributions 
to  viscous  shear  stress  at  the  wall.  Thus,  in  general,  a  perturbation 
analysis  of  a  boundary  layer  will  be  improper  near  a  point  where  the 
shear  stress  of  the  basic  flow  vanishes,  that  is,  near  a  separation 
point . 


Solution  of  Perturbation  Equations 

The  boundary- value  problem  for  fQ  (equations  (18)  and  (21))  arises 

in  the  case  of  boundary- layer  flow  on  a  flat  plate  and  has  been  studied 
thoroughly  in  the  literature .  For  completeness,  the  values  of  fQ  and 

its  derivatives,  as  determined  by  calculations  made  at  the  Lewis 
laboratory,  are  presented  in  table  I.  The  results  agree  precisely  with 
those  presented  in  reference  9 . 

The  functions  fj(X)  and  g-^(X)  may  be  expressed  as  linear  com¬ 
binations  of  functions  not  depending  on  the  parameters  M  or  Ag/fiA-^ . 

The  following  representation  was  proposed  by  L.  Richard  Turner  of  the 
Lewis  laboratory: 

gqCx)  =  f0(X)  +(l  +  Ili  M2)hx(X) 


(27a) 
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V 


(27b) 

The  coefficients  in  equations  (27)  may  he  determined  from  the  tables  of 
reference  3.  Substituting  equations  (27)  into  equations  (19),  (22), 
(23),  and  (25)  and  mating  use  of  equations  (18)  and  (2l)  yield  the 
following  differential  equations  and  boundary  conditions  for  the  func¬ 
tions  hx,  h2,  and  h3,  which  can  be  solved  once  and  for  all,  with  no 

regard  for  the  physical  parameters  of  the  problem: 


3*0*1"  -  2f0'h1'  +  eh-L- 


(28a) 


hi'(-)  =  h1'  (0)  =  h^O)  =  0  (28b) 

3f0h2"  +  3f0'’h2  +  6h2*"  =  0  (29a) 

h2'(»)  =1  h2'(0)  =  h2(0)  =  0  (29b) 

3f  h  "  +  3f  "h  +  6h  ' "  =  f  "h  (30a) 

0  3  0  3  3  0  1 


^•(-J  =  h3'(0)  =  h3(0)  =  0 


(30b) 


In  appendix  B,  the  computation  of  h-p  hg^  and  h^  is  described.  The 
results  appear  in  table  I. 


* 


Presentation  and  Discussion  of  Results 

The  various  important  boundary- layer  quantities  may  be  found  from 
the  profile  functions  given  in  equations  (27).  The  numerical  examples 
presented  in  figures  3  to  15  have  been  computed  with  the  aid  of  the 
tables  of  references  2  and  3. 

The  boundary  layer  on  a  cone  at  angle  of  attack  is  characterized 
by  the  phenomenon  of  "secondary  flow, "  which  has  been  described  by 
various  authors.  A  discussion  of  secondary  flow  as  it  occurs  on  a  cone 
at  angle  of  attack  is  provided  in  the  following  paragraphs  in  order  to 
facilitate  interpretation  of  the  results  of  the  present  analysis. 
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Since  the  inviscid  flow  about  the  cone  at  angle  of  attach  is 
"conical^ "  the  pressure  gradient  at  the  outer  edge  of  the  boundary  layer 
is  entirely  circumferential  (equations  (2)) .  The  streamlines  in  the 
inviscid  flow  just  outside  the  boundary  layer  (shown  schematically  as 
solid  lines  in  fig.  2)  will  incline  toward  the  direction  of  the  pres¬ 
sure  gradient .  Since  the  pressure  tends  to  be  constant  across  the 
boundary  layer;  the  fluid  in  the  boundary  layer  is  subject  to  the  same 
pressure  gradient  as  is  the  outer  flow;  but  has  less  inertia  with  which 
to  resist  its  effect;  and  thus  it  tends  to  follow  the  direction  of  the 
pressure  gradient  more  closely  than  does  the  outer  flow. 


An  angular  divergence  therefore  exists  between  outer  streamlines 
and  boundary-layer  streamlines;  the  divergence  is  greatest  when  the 
limiting  streamlines  at  the  base  of  the  boundary  layer  are  considered. 
These  limiting  streamlines  are  shown  schematically  as  dashed  lines  in 
figure  2.  The  divergence  e  is  defined  as  the  local  angle  included 
between  the  outer  and  limiting  streamlines . 

This  divergence  between  inner  and  outer  streamlines  ("secondary 
flow")  clearly  leads  to  a  draining  away  of  low-energy  air  from  the 
high-pressure  side  of  the  body  and  a  corresponding  concentration  of 
low-energy  air  on  the  low-pressure  side.  Accordingly;  if  the  cone  is 
at  positive  angle  of  attack;  the  skin  friction  would  be  expected  to  be 
lower  at  the  top  of  the  cone  than  at  the  bottom;  and  the  boundary- layer 
thickness  would  have  a  maximum  at  the  top  and  a  minimum  at  the  bottom. 


Before  proceeding  further;  it  is  perhaps  well  to  emphasize  that  the 
present  analysis  is  valid  only  in  the  limit  of  vanishing  angle  of  attack. 
That  is;  terms  multiplied  by  a  represent  rates  of  change  with  angle  of 
attack;  evaluated  at  zero  angle  of  attack.  All  results  will  be  presented 
in  this  form.  Whether  or  not  the  absolute  change  for  a  small  finite 
angle  of  attack  can  be  obtained  from  this  theory  depends  on  the  relative 
magnitude  of  the  effect  so  computed;  rather  than  on  the  size  of  a,  and 
depends  further  on  the  (unknown)  second  and  higher  derivatives  of  the 
quantity  with  respect  to  a. 


Velocity  profiles  and  streamline  divergence.  -  From  equations  (5} 
and  (leTH 


w  =  g\  =  a  A2  sin  <p  g± ' 

► 

u  =  fy  =  fQ*  -  a  Ax  cos  <p  f  » 


(31) 


These  expressions  may  be  divided  by  the  values  of  w-|_  and  ux  given 
in  equations  (2)  to  yield;  for  a«l ; 


w 


=  Si* 


(32a) 
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-^  =  f0'  +  a  Ax  cos  <p  (f0’  -  f^) 


(32b) 


or,  substituting  equations  (27), 


^  =  V  +11  +1LTMJhi 


(33a) 


£  -  v +  “  “1 005  ’ll1  - 1  s§(v  -  v>  -2!1  +  r-T  r)^  h5' 


-1  rA  A2 


(33b) 


The  outer  streamline  is  inclined  to  the  x-direction  by  an  angle 


-1  W1  .  +Bn-1  a  A2  sin  * 


tan_J-  =  tan 


1  -  a  cos 


a  Ag  sin  q> 


which  is  obtained  by  use  of  equations  (2)  .  The  limiting  streamline 
is  inclined  to  the  x-direction  by  an  angle 


tan_1[lim  zVa  iL  sin  <p 

\v*o  ^  f0”(o) 


which  is  obtained  with  the  aid  of  equations  (3l)  and  LfHospitalls 
rule.  The  divergence  6  is  given  by  the  difference  between  these  two 
angles: 


e  *  a  Ag  sin  cp 


8/(0) 

f0"(o) 


or,  by  use  of  equation  (27a)  and  table  I, 


c  «  a(l.29l)  Ag^L  +  M^sin  <P 


This  streamline  deflection  due  to  angle  of  attack  is  indicated  by  the 
quantity 
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— r— —  =  (1.291)  A2  fl  +  2li  M2) 

am  <P  W/^o  V  2  1 

and  is  shown  in  figure  3  as  a  function  of  stream  Mach  number  Mq  for 
cone  semivertex  angles  ©  of  10°,  20°,  and  30°. 

The  distribution  through  the  boundary  layer  of  the  circumferential 
velocity  w/w-^  is  shown  in  figure  4  for  cone  semivertex  angles  ©  of 

10°  and  30°,  and  for  stream  Mach  numbers  Mq  of  2  and  4.  From  equa¬ 
tion  (31)  and  figure  4,  it  is  evident  that  no  "cross-flow  separation" 
takes  place,  because  the  profiles  of  figure  4  (which  are  not  of  the 
separation  type)  appear  at  all  circumferential  locations,  modified  by 
the  value  of  sin  <P,  which  affects  only  the  scale  of  the  w-profiles . 

Of  course,  cross-flow  separation  would  occur  for  large  angles  of 
attack  because  of  the  effect  of  an  adverse  circumferential  pressure 
gradient  over  part  of  the  cone  surface.  For  very  small  angles  of 
attack,  separation  of  the  cross  flow  would  hardly  be  expected  since  the 
circumferential  pressure  gradient  is  always  favorable  (equations  (2)). 

The  increment  in  meridional  velocity  profile  due  to  angle  of  attack 
is  indicated  by  the  quantity 

fe]L0  =  \  (V-V} 

which  is  shown  in  figure  5  for  ©  =  10°  and  30°,  and  for  -M0  =  2  and 
4.  For  positive  angle  of  attack,  the  meridional  profile  is  correspond¬ 
ingly  steeper  on  the  underside  and  less  steep  on  the  top.  This  result 
agrees  with  the  qualitative  argument  that  low-energy  air  should  drain 
away  from  beneath  the  body  and  concentrate  on  the  top.  The  order  of 
magnitude  of  this  effect  is  shown  in  figure  6,  in  which  the  u/u-j^  pro¬ 
files  in  the  plane  <p  =  0,Jt  for  a  =  2.5°  are  compared  with  the  pro¬ 
files  for  a  =  0,  when  ©  =  10°  and  Mq  =  4. 

The  relative  thickening  of  the  boundary  layer  at  the  top  of  the 
cone  and  thinning  at  the  bottom  might  be  expected  to'  have  an  adverse 
effect  on  laminar  stability  at  the  top  and  the  reverse  effect  on  the 
bottom . 

Skin  friction.  -  The  circumferential  and  meridional  components  of 
the  viscous  shear  stress  at  the  cone  surface  may  be  written  as 
coefficients: 
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A I  l  A-2 

C~  W—  =  0.664  1  -  a  A,  cos  q>  1.502  -  0.525  - 

rxy  3C  X  V  0A-L 


y-X  —2  ^2 

0.190  —  M  - 


2A^ 


(35b) 


The  meridional  skin  friction  at  zero  angle  of  attack  is  obtained 
by  setting  a  =  0  in  equation  (35b),  yielding  a  result  equivalent  to 
the  results  of  Hantsche  and  Wendt  (reference  6)  and  Mangier. 

In  figure  7,  the  circumferential  skin  friction  appears  as  a  func¬ 
tion  of  stream  Mach  number  for  semi vertex  angles  of  10°,  20°,  and  30° 


Figure  8  shows  the  increment  due  to  angle  of  attack  of  meridional 
skin  friction  plotted  as  a  function  of  Mo  for  the  same  cone  angles. 

The  effect  of  angle  of  attack  on  meridional  skin  friction  is  of  the  same 
order  of  magnitude  as  the  effect  on  circumferential  skin  friction,  as 
would  be  expected  in  view  of  the  previous  discussion  of  the  velocity  pro¬ 
files.  For  positive  angle  of  attack,  the  meridional  skin  friction  has  a 
maximum  at  <p  =  0  and  a  minimum  at  q>  =  it,  which  is  consistent  with  the 
previous  qualitative-discussion  of  secondary  flow. 

Boundary- layer  thickness .  -  The  boundary- layer  thicknesses  are 
expressed  as  the  mass-flow  defects  ("displacement  thicknesses") 
associated  with  the  u  and  w  profiles .  Thus  i 


6<P  - 


1  -  Pv 
p  w 

i  1  1 


or,  making  use  of  equations  (l),  (7),  and  (8)  and  expressing  8  as  a 
Reynolds  number  R$  =  p  u  S/p, 
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1/  3 


x/P\  P 


x  -#e=)2 i 


h  -  A 

\P  vj 


Pi  _  u_ 
P  '  U, 


From  equations  (9c)  and  (9d), 


p  p/p  2  Lu2 


Substituting  equations  (5)  and  (37)  into  (36), 


’«p  II  3  W  Pi 


s  £-  1+  ni  if  .  h! 

O/h  I  %  \  TTk  I  1T\  / -Y\  .  O  9  I  4 


p/p  2  _2 

'  u 


2/  w- 


-  —  dX  (38a) 


’x  U  3 


H. .  =raf  ?-  1  ♦  P&  ^  S2  ^-[1  -  4V  i  ax  (38b) 


p/p  2  u2  V  U  2/  U1 


Introducing  equations  (2)  and  (33a) ,  performing  the  indicated  integrations,, 
and  making  use  of  table  I  change  equation1  (38a)  to 


R-.  A  -=r-  =  0.944  +  1.608  m2 
°<P  WCRX  d 


(39a) 


It  would  be  improper  to  present  the  term  of  order  a  resulting  from 
integration  of  equation  (38a)  since  in  the  present  approximation  the 
ratio  w/w1  is  known  only  to  unit  order  (see  equation  (33a)).  ,  Intro¬ 
ducing  equations  (2)  and  (33b)  into  equation  (38b)  and  performing  the 
required  integrations  (with  the  aid  of  equations  (18),  (19),  and  (25)) 
result  in  the  following  expression: 
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=  1.721  +  2.385  Eli  M2 


1  +  a  A,  cos  <p 


1.721  +  2.385  M2 


0.863  - 


3.576  M2 


-  <0.863  + 


1.469  Eli  m2  + 


rr_i  _2  2\  a2 

0.385  V^M  >  — 

L  2  J  J  SAj^ 


2.385 


-1  r2)h 


M  S—  - 


(39b) 


The  displacement  thickness  associated  with  the  u-conrponent  at  zero 
angle  of  attack  (obtained  by  setting  a  =  0  in  equation  (39b))  is  shown 
in  figure  9,  as  a  function  of  for  various  cone  angles,  and  is 

equivalent  to  the  results  of  Hantsche  and  Wendt  (reference  6)  and  Mangier. 

The  w-thickness  given  by  equation  (39a)  is  shown  in  figure  10.  This 
thickness  does  not  depend  on  <p,  and  hence,  at  any  value  of  x,  is  con¬ 
stant  around  the  cone. 

The  increment  in  u-thickness  due  to  angle  of  attack  (fig.  ll) 
indicates  thinning  of  the  boundary  layer  on  the  high-pressure  side  of 
the  cone  and  thickening  on  the  low-pressure  side.  This  effect  is  qual¬ 
itatively  consistent  with  the  previous  discussion  of  secondary  flow  in 
the  boundary  layer  of  a  cone  at  angle  of  attack,  and  its  order  of 
magnitude  is  consistent  with  the  corresponding  results  for  velocity 
profiles  and  skin  friction. 

Forces  and  moments.  -  The  inviscid  solution  for  flow  about  the  cone 
at  an  angle  of  attack  provides  the  forces  and  moments  due  to  the  action 
of  normal  pressure  forces,  which  may  be  evaluated  using  the  tables  of 
references  2  and  3.  In  order  to  obtain  the  total  force  and  moments,  the 
forces  and  moments  due  to  viscous  shear  at  the  surface  must  be  added  to 
the  corresponding  pressure  forces  and  moments.  Lift  and  drag  are  measured 
normal  and  parallel,  respectively,  to  the  stream  direction.  Pitching 
moment  is  measured  about  the  cone  apex,  positive  in  the  sense  shown 
in  figure  12 .  All  quantities  appearing  in  the  following  formulas  are  in 
dimensional  form. 
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(a)  Viscous  drag:  Written  in  coefficient  form,  the  drag  due  to 
skin  friction  at  the.  surface  is 


C  — 


v  1 


Dv 

— 2 

1  pu 

-  p  u  2  A 

2  H0  0 

A  Po-V 

ix  p2« 


dx  cf  eyi  -  e*  x  d«p 

X 


where  D  is  the  viscous  drag  and  A  is  the  cone  surface  area.  Using 
equation7 (55b)  yields 


C  =  0.885/1  P  —  —  ^  -  e2 


Rx  p0  Mq2 


This  expression  gives  the  drag  at  zero  angle  of  attack.  The  correction 
due  to  angle  of  attack  is  of  order  a2 ,  and  hence  is  neglected  in  the 
present  linear  analysis.  The  drag  coefficient  (equation  (40.)),  which 
is  equivalent  to  the  results  obtained  by  Hantsche  and  Wendt  (reference  6) 
and  Mangier,  is  shown  in  figure  13  as  a  function  of  Mq  for  various 

values  of  ©  * 

(b)  Viscous  lift:  The  lift  due  to  circumferential  skin  friction 


C  s  - - - 

Lv  -1  2 

ip0u0  A 


=  -a  CD  +  -  -2^ 
Dv  A  5 


dx  0x  sincp  d<p 


where  is  the  viscous  lift.  Using  equation  (35a), 

CT  =  a  -  C-  +  (l.015  +  0.572  Ui  M^A-  E__-lOp 

\  Dv  V  2  7  2  P0  ^  K 


The  lift  coefficient  is  shown  in  figure  14  as  a  function  of  Mq  for 
various  values  of  ©  .  Except  for  quite  low  Reynolds  numbers,  the  laminar 
viscous  lift  is  small  compared  with  the  pressure  lift.  The  pressure 
lift  coefficient  is  of  order  a,  whereas  for  a  Reynolds  number  based  on 
slant  height  of  300,000,  the  viscous  lift  coefficient  on  a  10°  cone  at 
Mq  =  4  is  only  of  the  order  of  0.01  a; 
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(c)  Viscous  pitching  moment:  By  use  of  equation  (35a),  the  moment 
due  to  circumferential  skin  friction  may  he  found : 


(42) 

where  m^  is  the  viscous  pitching  moment  about  the  apex.  The  moment 
is  shown  in  figure  15  as  a  function  of  Mq  for  several  cone  angles. 


CONCLUSIONS 

For  small  angles  of  attack,  the  outer  nonviscous  flow  about  a  cone 
at  angle  of  attack  to  a  supersonic  stream  has  previously  been  treated 
by  a  perturbation  theory,  which  may  be  used  to  describe  the  outer  flow 
governing  the  development  of  the  laminar  boundary  layer  on  the  cone 
surface,  provided  that  no  account  need  be  taken  of  the  "vortical  layer" 
recently  discussed  by  Ferri. 

This  vortical  layer  is  a  region  in  the  nonviscous  flow  field  near 
the  surface  of  the  cone  across  which  there  are  large  gradients  of 
entropy  and  velocity,  and  is  said  to  have  an  angular  thickness  of  the 
order  of  the  square  of  the  angle  of  attack.  If  the  boundary  layer  is 
much  thicker  than  this,  then  presumably  viscous  diffusion  disperses  the 
vortical  layer,  which  then  may  be  neglected. 

Since  the  outer  nonviscous  flow  is  "conical, "  use  may  be  made  of 
previously  derived  differential  equations  for  a  boundary  layer  governed 
by  an  outer  flow  with  conical  symmetry.  These  equations  imply  a  solution 
having  parabolic  similarity  in  meridional  planes .  Assuming  a  Prandtl 
number  of  1,  no  heat  transfer  through  the  surface,  and  validity  of  a 
linear  temperature-viscosity  relation,  it  is  necessary  only  to  solve  the 
differential  equations  for  two  components  of  the  vector  potential. 

In  the  limit  of  vanishing  cone  angle  (finite  angle  of  attack),  the 
boundary-layer  equations  become  precisely  those  governing  the  compressible 
flow  over  a  yawed  infinite  cylinder  of  diameter  equal  to  the  local  diameter 
tke  cone  at  the  point  under  consideration.  The  cross  and  spanwise 
flows  are  coupled  through  the  density,  however,  and  it  is  not  proper  to 
compute  the  cross  flow  by  two-dimensional  methods,  as  may  be  done  in  the 
incompressible  case.  This  coupling  depends  on  the  stream  Mach  number 
and  not  on  the  Mach  number  of  the  cross  flow. 
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In  the  limit  of  vanishing  angle  of  attack  (finite  cone  angle)  ,  the 
■boundary- layer  equations  reduce  to  those  governing  flow  about  a  cone  at  * 

zero  angle  of  attack.  A  sufficiently  small  angle  of  attack  may  he  con¬ 
sidered  to  impose  small  perturbations  on  the  flow  at  zero  angle  of  attack. 

If  these  perturbations  are  written  in  a  form  appropriate  to  the  stream 
boundary  conditions  provided  by  the  nonviscous  perturbation  theory, 
ordinary  differential  equations  result.  These  may  then  be  solved  by 
numerical  integration  to  obtain  the  fluid  velocities  and  state  quantities 
in  the  boundary  layer.  The  type  of  perturbation  method  used  in  this 
report  would  be  improper  in  the  vicinity  of  any  separation  point  of  the 
basic  flow.  No  such  difficulty  arises  in  the  present  application. 

According  "to  the  perturbation  analysis,  the  boundary-layer  thick¬ 
ness  associated  with  the  cross  (circumferential)  velocity  is  constant 
around  the  body  at  a  given  distance  from  the  apex.  The  shear  associated 
with  this  cross  flow  gives  rise  to  lift  and  moment  on  the  body.  No 
separation  of  the  cross  flow  occurs,  presumably  as  a  consequence  of  the 
favorable  circumferential  pressure  gradient. 

Secondary  flow  occurs  in  the  boundary  layer  of  a  cone  at  angle  of 
attack;  the  slow  fluid  near  the  base  of  the  boundary  layer  tends  to 
incline  more  toward  the  direction  of  the  circumferential  pressure 
gradient  then*  does  the  outer  flow.  This  mechanism  causes  a  concentration 
of  low-energy  air  at  the  low-pressure  side  of  the  cone  (top,  for  positive 
angle  of  attack)  .  Therefore,  taking  into  account  the  change  in  the  pro¬ 
file  of  meridional  velocity  due  to  angle  of  attack,  the  boundary-layer 
thickness  associated  with  this  profile  is  greater  at  the  top  and  less  on 
the  bottom  of  the  cone  than  the  value  appropriate  to  the  case  of  zero 
angle  of  attack.  The  reverse  is  true  of  meridional  skin  friction. 

Meridional  skin  friction  is  integrated  to  give  the  viscous  drag  of  the 
cone.  The  effect  of  angle  of  attack  on  drag  is  second  order  in  angle 
of  attack  and  is  ignored  in  the  present  analysis. 

The  changes  in  the  meridional  velocity  profile  due  to  angle  of 
attack  suggest  that  a  decrease  in  laminar  stability  on  the  low-pressure 
side  and  the  reverse  effect  on  the  high-pressure  side  might  be  expected. 

Lewis  Flight  Propulsion  Laboratory 

National  Advisory  Committee  for  Aeronautics 
Cleveland,  Ohio,  August  10,  1951 
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APPENDIX  A 

SYMBOLS 


*1 

AZ 

\ 


Cmv 

CP 

Dv 

f 

fO> 

8 

gl 


The  following  symbols  are  used  in  this  report: 

related  to  perturbation  in  meridional  velocity  u  outside 
boundary  layer 

related  to  perturbation  in  circumferential  velocity  outside 
boundary  layer 

related  to  perturbation  in  pressure  outside  boundary  layer 

related  to  perturbation  in  density  outside  boundary  layer 

constant  of  proportionality  in  temperature-viscosity 
relation 

coefficient  of  viscous  drag 
coefficient  of  viscous  lift 

coefficient  of  viscous  pitching  moment 
specific  heat  at  constant  pressure 
viscous  drag 

component  of  vector  potential  related  to  meridional 
velocity  u 

f]_  value  of  f  for  a  =  0  and  perturbation  of  f , 

respectively 

component  of  vector  potential  related  to  circumferential 
velocity 

linear  approximation  to  g 
h0,  bu  functions  composing  fn  and  gn 
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Subscripts : 

w  denotes  evaluation  at  cone  surface 

0  denotes  evaluation  in  undisturbed  stream 

1  denotes  evaluation  at  outer  edge  of  boundary  layer 

Subscript  notation  for  partial  differentiation  is  used. 

Superscripts : 

1  Primes  denote  ordinary  differentiation, 

~  Bars  (as  in  u  or  R)  denote  evaluation  in  nonviscous 

flow  next  to  cone  surface  in* case  of  zero  angle  of 
attack , 
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APPENDIX  B 

NUMERICAL  DETERMINATION  OF  PROFILE  FUNCTIONS 
By  L .  Richard  Turner 

The  four  functions  fQ(X),  h^X),  hg(X),  and  h3(X)  have  "been 

obtained  by  numerical  solution  of  the  boundary- value  problems  set  forth 
in  equations  (18)  and  (28)  to  (30) . 

Differential  equation  (18)  is  the  basic  Blasius  equation*  it  has 
been  solved  to  provide  the  values  of  the  function  fQ(X)  required  for 

the  remaining  integrations.  Boundary  conditions  (2l)  are  specified  at 
two  points  X  =  0 ,<a  and  thus  cannot  be  satisfied  directly  by  a  step- 
by-step  integration  process.  This  difficulty  may  be  avoided  as  follows: 
Define 

q(a)  =  afQ(\) 

a  s  -  \ 
a 


■where  a  is  an  arbitrary  constant.  Substitution  into  equations  (18) 
and  (2l)  yields 

q  q"  +  2q ,n  =  0  (Bl) 

ql(«>)  -  a2  q’(0)  =  q(0)  =  0  (B2) 

Differential  equation  (Bl)  may  he  integrated  using  the  boundary  con¬ 
ditions  q’(0)  =  q(o)  =  0  and  an  arbitrary  value  of  q"(o).  A  constant 
value  of  q*  (a)  will  ultimately  be  reached  and  may  be  equated  to  a2 
since  a2  is  arbitrary.  The  resulting  value  of  a  may  then  be  used  to 
invert  the  transformation  leading  to  equations  (Bl)  and  (B2) ,  thus 
determining  the  solution  Tq(X)  • 

Equations  (28)  to  (30)  are  linear  in  their  dependent  variables; 
solutions  that  satisfy  the  boundary  conditions  at  X  =  *>  may  therefore 
be  found  by  appropriate  linear  combinations  of  the  complementary  and 
particular  solutions  satisfying  the  given  boundary  conditions  at  the 
origin . 

The  functions  h^,  h2,  and  h3  were  found  by  computing  the 

complementary  and  particular  integrals  of  equations  (28a)  and  (30a) 
subject  to  the  following  boundary  conditions: 
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Equation 

Integral 

h(0),  h'(0) 

h"(0) 

(28a) 

Complement  ary 

0 

1/3 

(28a) 

Particular 

0 

1/3 

(30a) 

Complementary 

0 

1/3 

(30a) 

Particular 

0 

0 

The  values  of  h"(o)  are  arbitrary. 

The  function  h-j_  is  thus  composed  of  the  complementary  and 
particular  solutions  of  equation  (28a),  linearly  combined  in  such  a 
way  that  h-j/ (•)  =  0  (boundary  condition  (28b)).  The  function  h2(X) 

was  obtained  by  scaling  the  complementary  integral  of  equation  (30a), 
■which  satisfies  equation  (29a),  so  that  h^*)  =  1  (boundary 

condition  (29b)).  The  function  hj(X)  is  composed  of  the  complementary 
and  particular  integrals  of  equation  (30a),  combined  linearly  so  that 
h^ 1  (•)  =  0  (boundary  condition  (30b)). 

The  four  differential  equations  (18a)  and  (28a)  to  (30a)  were 
integrated  numerically  by  a  method  somewhat  similar  to  that  of  Adams 
(see  reference  ll) .  Given  the  complete  solution  for  the  function  in 
question  at  six  equally  spaced  points,  the  solution  may  be  extended  to 
the  next  point  as  follows:  A  sixth-degree  polynomial  is  passed  through 
the  six  known  and  one  unknown  values  of  the  third  derivative.  This 
polynomial  may  be  integrated  to  give  the  value  and  first  two  derivatives 
of  the  function  at  the  seventh  point  in  terms  of  the  unknown  third 
derivative.  Substitution  of  these  results  into  the  differential  equation 
evaluated  at  the  seventh  point  yields  an  algebraic  equation  for  the 
third  derivative  at  that  point,  which  may  be  solved,  thus  extending  the 
solution  of  the  differential  equation  by  one  step.  This  process  is 
repeated  out  to  sufficiently  large  values  of  X  so  that  the  boundary 
conditions  at  X  =  «*  may  be  satisfied  to  the  desired  degree  of  accuracy. 
A  step  size  of  0.01  was  used  for  equation  (18a),  and  a  step  size  of  0.02 
was  used  for  equations  (28a)  to  (30a). 


The  use  of  this  method  requires  that  an  initial  set  of  values  of 
each  function  and  its  first  two  derivatives  be  known  in  advance  for  the 
first  six  intervals,  starting  at  the  origin.  These  values  were  obtained 
by  use  of  a  method  of  successive  approximation  which  is  a  numerical 
application  of  the  classical  Picard  iteration  process. 

All  the  numerical  integrations  were  performed  on  the  IBM  Card 
Programmed  Electronic  Calculator,  and  carried  out  to  ten  decimal  places. 
The  computations  were  checked  by  computing  the  third  and  seventh 
differences  of  the  third  derivatives  to  assure  freedom  from  computational 
error  and  from  significant  truncation  error,  respectively. 
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VO 

CO 

00 

00 


<\ 

-0.031614 

-.031570 

-.031270 

-.030488 

-.029042 

-.026801 

-.023700 

-.019750 

-.015043 

-.009752 

-.004120 

.001563 

.006985 

.011845 

.015887 

.018928 

.020872 

.021720 

.021559 

.020544 

.018877 

.016775 

.014450 

.012086 

.009829 

.007780 

.006000 

.004511 

.003308 
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.000755 

.000493 
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.000013 
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.  000004; 
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.OOOOOlj 

.000001 

0 

w 

to 

A 

) 

-.006321 

-.012615 

-.018810 

-.024794 

-.030421 

-.035525 

-.039934 

-.043485 

-.046039 

-.047502 

.047827 

.047034 

.045201 

.042463 

.039002 

.035028 

.030762 

.026417 

.022194 

.018243 

.014673 

.011548 

.008895 

.006706 

.004949 

.003576 

.002530 

.001753 

.001189 

.000790 

.000514 

.000328 

.000205 

.000125 

.000075 

.000044 

.000025 

.000014 

.000008 

000004 

000002 

000001 

000001 

-p- 

1  1  1  1  1 

t  1  t  1  1 

1  1  i  t  t 

<1111 

1  1  i  1  1 

*  1  1  1  0 

0  0 

J? 

'  0 

-.000632 

-.002526 

-.005671 

-.010035 

-.015562 

-.022165 

-.029721 

-.038075 

-.047041 

-.056410 

-.065957 

-.075457 

-.084693 

-.093470 

-.101625 

-.109033 

-.115614 

-.121331 

-.126189 

-.130227 

-.133512 

-.136126 

-.138163 

-.139715 

-.140874 

-.141721 

-.142326 

-.142750' 
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Figure  5.  -  Increment  in  meridional  profile  due  to  angle  of  attack. 
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Figure  7.  -  Circumferential  skin  friction  due  to  angle  of  attack. 
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Figure  11.  -  Increment  in  displacement  thickness  of  meridional  velocity- 
profile  due  to  angle  of  attack. 
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Figure  15.  -  Coefficient  of  viscous  pitching  moment  about  done  apex. 
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